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We theoretically study a competition between two types of spin cycloids (a6-plane and de- 
plane ones) in the multiferroic perovskite manganites, which is an origin of intriguing magne- 
toelectric phenomena in these compounds. Analysis of a microscopic model using the Monte- 
Carlo method reveals that their competition originates from a conflict between the single-ion 
anisotropy and the Dzyaloshinsky-Moriya interaction. It is demonstrated that the conflict can 
be controlled by tuning the second-neighbor spin exchanges through the GdFeOs-type distor- 
tion, which leads to a cycloidal-plane flop from ab to be observed in the solid-solution systems 
like Eui-xYccMnOs and Gdi-xTbccMnOs with increasing x. 
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Magnetoelectric (ME) coupling in concurrently mag- 
netic and ferroelectric [multiferroic (MP)] materials 
causes a lot of fascinating ME phenomena.^' ^ Inten- 
sive experiments revealed magnetic-field {H) induced 
ferroelectric-polarization (P) flops, ^ giant magnetoca- 
pacitance effect s,^^ and electrically activated magnon 
excitations (termed electromagnons)^'^ in the typical 
MF materials of perovskite manganites i^MnOa with R 
being a rare-earth ion. The coupling and the phenom- 
ena have attracted enormous interest from viewpoints of 
both science and application since their manipulations 
enable us magnetic control of ferroelectricity and elec- 
tric control of magnetism. 

The ferroelectricity in i^MnOa is induced by the 
cycloidally ordered spin structure through the inverse 
Dzyaloshinsky-Moriya (DM) mechanism. The de- 
plane (a6-plane) spin cycloid propagating along h gen- 
erates P along the c (a) axis.^^ Since the first report 
of MF nature in TbMnOs, a number of MF materials 
of spiral-magnetism origin have been discovered in suc- 
cession, and these materials have turned out to exhibit 
similar interesting ME phenomena also. However, up to 
now, mechanisms of these ME phenomena have not been 
clarified yet in spite of much effort. 

An important clue to approach these issues were re- 
cently provided from precise experimental studies on the 
phase diagrams of i^MnOs. When we vary the mag- 
nitude of GdFeOa-type distortion continuously by us- 
ing solid solutions like Eui-^Y^MnOa (EYMO)^^'^^ and 
Gdi-^Tb^^MnOs (GTMO),^"^ the 90° reorientation of P 
from P||a to P||c is observed as Y or Tb concentration x 
increases, which is accompanied by a cycloidal-plane flop 
from ah to hc^^ see Fig. 1. 

Interestingly we notice that the competition between 
these two types of spin cycloids {ah /he competition) has 
a universal relevance to the above ME phenomena. The 
P flop under an applied is a consequence of the compe- 
tition controlled by H. The theoretically proposed elec- 
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Fig. 1. (Color online) Experimental T-x phase diagrams of (a) 
Eui-xYccMnOs and (b) Gdi-^TbccMnOs reproduced from 
Ref.-*^^ and Ref./"^ respectively. PE and FE denote paraelectric 
and ferroelectric phases, and AFM(A)+WFM denotes the A- 
type antiferromagnetic phase with weak ferromagnetism. 



tromagnon is an oscillation of the spin-cycloidal plane. 
The giant magnetocapacitance effect was experimentally 
ascribed to the electric-field-driven motion of domain 
walls between the ah- and 6c-cycloidal phases.^ 

These indicate that clarification of origin of the ah /he 
competition in i^MnOs is a key to understanding the ME 
phenomena in these MF materials. Recent theoretical 
studies reproduced the spiral spin order in i^MnOa.^'^^ 
However the origin of the ah /he competition and a mech- 
anism of the cycloidal-plane flop were not addressed. 

In this Letter, by taking i^MnOs as a typical exam- 
ple, we theoretically study the ah /he competition by con- 
structing a microscopic model for the Mn Sd-spin sys- 
tem. We reveal that the competition originates from a 
conflict between the single-ion anisotropy and the DM 
interaction. This competition is controlled by the second- 
neighbor spin exchanges enhanced by the GdFeOs-type 
distortion, which leads to a cycloidal-plane flop in EYMO 
and GTMO with increasing x. We demonstrate that their 
T-x phase diagrams are naturally reproduced as its con- 
sequence. Our finding provides a fundamental insight on 
the ME phenomena in the MF materials of magnetic- 
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Fig. 2. (Color online) Superexchange interactions, tilted local 
axes ^i, r]i, and attached to the i-th MnOe octahedron, and 
Dzyaloshinsky-Moriya vectors in i^MnOs. 

spiral origins. 

We start with a classical Heisenberg model with some 
additional interactions and magnetic anisotropics on a 
cubic lattice, in which the Mn S=2 spins are treated as 
classical vectors. The Hamiltonian consists of four terms 
as = Hex + '^^sia + '^^DM + '^^cub-^^ We adopt meV as 
the unit of energy. 

The first term TYex denotes superexchange interac- 
tions; ferromagnetic (FM) Jab on the in-plane nearest- 
neighbor bonds, antiferromagnetic (AFM) J2 on the in- 
plane second-neighbor bonds along the b axis, and AFM 
Jc on the bonds along the c axis as shown in Fig. 2. The 
second-neighbor exchanges J2 are caused by finite over- 
laps of Cg orbitals in the b direction mediated by two 
oxygens. 

The second term Tigia denotes the single-ion anisotropy 
(SIA), which can be written as TYgia = ^'^.{Si ■ 

Here ^i, r]i and d are directional vectors of the tilted 
local axes attached to i-th MnOe octahedron as shown 
in Fig. 2. Because of this term, the c axis in i^MnOa 
becomes a hard magnetization axis. For the vectors, we 
use the structural data of EuMnOa.^^ We have confirmed 
that the results are not significantly changed even if we 
use structure data of other i^MnOa. 

We have microscopically determined all the parame- 
ters in this model. The spin exchanges Jab and Jc and 
the SIA parameters D and E have been calculated by us- 
ing formulae given in Ref.^^ and Ref.,^^ respectively. We 
have found that the values are nearly independent of R as 
far as vicinities of the multiferroic phases are concerned. 
We take 7^6=0.80, Jc=1.25, D=0.25, and E=0.30. 

We study the x dependence by varying the value of 
J2. We consider that main roles of the GdFeOs-type dis- 
tortion controlled by x on the magnetic properties in 
i^MnOs are induction and enhancement of the second- 
neighbor antiferromagnetic exchanges J2.^^ This is be- 
cause their exchange paths contain two oxygen 2p or- 
bitals, and the distortion enhances the hybridization be- 
tween these two 2p orbitals. Consequently, the spiral ro- 
tation angle ^=180°xg'Mn with qun being a spiral wave 
number increases as x increases (note that a relation 
coscj) = Ji/(2J2) holds in the simple J1-J2 model). The 



crucial importance of the second-neighbor exchanges J2 
is peculiar in these Mn^+ compounds. Because of the 
tlg^l electron configuration, the nearest-neighbor FM 
coupling Jab is strongly reduced due to the cancellation 
of opposite contributions from two different orbital sec- 
tors, i.e. FM coupling between 5=1/2 Cg spins and AFM 
coupling between S=3/2 t2g spins. This leads to the in- 
creasing importance of the exchanges J2 with cooperative 
AFM contributions from these two orbital sectors. 

The third term Hbm = dij - {Si x Sj) denotes 

the DM interactions.^^' The DM vectors dij are de- 
fined on the Mn(i)-0-Mn(j) bonds, and are expressed in 
terms of five DM parameters, a^b, /3ab^ labi <^c and 
as shown in Fig. 2 reflecting the crystal symmetry. We 
can deduce their values from results of the first-principles 
calculation^^ and the electron-spin resonance (ESR) ex- 
periments.^^ We take aa6=0.10, /3a6=0.10, 7^6=0.14, 
Q/^=0.30, and ;^c=0.30. We also study the case for a 
slightly larger ac of 0.38 without changing the other pa- 
rameters. The last term Hcuh = s{s+i) ^i^^ti + + 
S^^) represents the cubic anisotropy. The constant a for 
Mn^+ ion in the Oe octahedron was evaluated to be 
0.0162 in the ESR experiment. 

We calculate thermodynamic properties of this model 
by using the Monte-Carlo (MC) method. To overcome 
the slow MC dynamics at low temperature (T), we em- 
ploy the replica exchange MC method. Each exchange 
sampling is taken after 400 standard MC steps. Typically, 
we perform 600 exchanges for a system with 48x48x6 
sites under the periodic boundary condition. We confirm 
that the finite-size effect is small enough and never affects 
our conclusion qualitatively. We calculate specific heat 
Cs{T) and total spin-helicity vector h^(T) to determine 
the transition points and the magnetic structures, which 
are respectively calculated as Cs{T) = j^d{l-L) / d{k^T)^ 
and h^{T) = j^{\Y.i^i^ Si^b\)/S'^, where the brack- 
ets denote thermal averages. The assignments of transi- 
tions and phases are confirmed by the calculated spin- 
correlation functions. 

In Figs. 3(a)-3(c), we show the calculated Cs{T) and 
h^{T) when ac=0.30 for various values of J2. Here h^{T) 
(7=a, 6, c) denotes the 7 component of h^{T). In the ab- 
cycloidal [6c-cycloidal] phase, hc{T) [ha{T)] has a large 
value, while other two components are strongly sup- 
pressed. On the other hand, in the AFM(A) and si- 
nusoidal collinear phases, all of the three components 
should be nearly equal to zero. 

In Fig. 3(a), we can see a single phase transition in 
Cs{T) for a small value of J2=0.36, which corresponds 
to a transition from paramagnetic to canted AFM (A) 
phases. Through this transition, h^{T) is nearly equal to 
zero constantly. 

For a larger value of J2=0.62, CsiT) shows two peaks 

indicative of two thermal transitions see Fig. 3(b). 

h^{T) is approximately zero constantly through the first 
transition, at which the system enters into the sinusoidal 
collinear phase from the paramagnetic phase. Contrast- 
ingly, at the subsequent transition into the a6-cycloidal 
phase, its c component hc{T) starts increasing. 

For a further increased value of J2=0.80, the sys- 
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Fig. 3. (Color online) Calculated temperature profiles of specific 
heat Cs{T) and total spin helicity h\T) = [ha(T), hb(T), hc(T)] 
for (a) J2=0.36, (b) J2=0.62, and (c) J2=0.80 when ac=0.30. 



tern undergoes three thermal transitions from high to 

low temperatures see Fig. 3(c). The first one is 

a transition from paramagnetic to sinusoidal collinear 
phases through which all of the three components of 
h^{T) are approximately zero constantly . At the sec- 
ond transition into the 6c-cycloidal phase, its a compo- 
nent ha{T) increases, while other two components re- 
main to be small. With further lowering T, the a compo- 
nent hair) suddenly drops, while the c component hc{T) 
steeply increases at the third transition accompanied by 
a cycloidal-plane flop from be to ab. 

In Figs. 4(a) and 4(b), we display theoretically ob- 
tained T-J2 diagrams for ac=0.30 and ac=0.38, respec- 
tively. They are in good agreement with the experimental 
T-x diagrams of EYMO and GTMO shown in Fig. 1. For 
both EYMO with x=OA and TbMnOs, the value of J2 
is calculated to be approximately 0.65 from experimen- 
tally obtained ^Mn-^^'^^ Both diagrams show that the 
6c-cycloidal regime increases with increasing J2, result- 
ing in the cycloidal-plane flop in certain T ranges. 

Here we argue the mechanism of the flop. In the ab- 
cycloidal state the rotating spins couple dominantly to 
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Fig. 4. (Color online) Theoretical T-J2 phase diagrams for (a) 
q;c=0.30 and (b) q;c=0.38. PE and FE denote paraelectric and 
ferroelectric phases expected in the inverse DM model, re- 
spectively. 



the c components of DM vectors on the in-plane Mn-0- 
Mn bonds. Their magnitudes are all equal to 7a6, and 
their signs (i.e. +7a6 and —^ah) are alternately arranged 

along the x and y bonds see Fig. 5(a). Without DM 

interaction, the spins rotate with the uniform rotation 
angles of (\)ah' On the other hand, in the presence of DM 
interactions, the rotation angles become to be alternately 
modulated into (\)ah^^<\>ah and (j)ab — ^4^ab with A(j)ab > 

to get an energy gain from the DM interactions see 

the inset of Fig. 5(a). We can derive the energy gain due 
to this angle modulation as 

AE^\^/N = --fabS^\sm{(t)ab - A(t)ab) - siu | 

= -7a6^^|cOS(/)a6|A(/)a6. (1) 

This expression implies that the energy gain |A£^gj^| 
is reduced with increasing (j)ab because the prefactor 
Icos^abl becomes maximum (=1) for (pab =0 but de- 
creases as (j)ab increases. Since the angle (j)ab increases 
as J2 increases, the a6-cycloidal state is destabilized with 
increasing J2 or with increasing GdFeOa-type distortion. 

On the other hand, the spins in the 6c-cycloidal state 
dominantly couple to the a components of DM vectors 
on the out-of-plane Mn-O-Mn bonds. Their magnitudes 
are all equal to ac, and their signs are the same within 

a plane, but alternate along the c axis see Fig. 5(b). 

Without DM interaction, the angles between adjacent 
two spins along the c axis are uniformly 0c = tt because 
of the strong AFM coupling Jc. In the presence of DM 
interactions, the angles again suffer from modulations 

into TT + A(j)c and tt — A(j)c with A^c > see the 

inset of Fig. 5(b). Similar to the a6-cycloidal case, we 
can derive the energy gain due to the angle modulation 
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Fig. 5. (Color online) (a) [(b)] Spin structure in the a6-cycloidal 
[6c-cycloidal] state and arrangement of the c [a] components of 
DM vectors on the in-plane [out-of-plane] Mn-O-Mn bonds. The 
symbols and (g) express their signs, i.e. for the positive sign 
and for the negative sign. Inset shows spin directions with 
modulated rotation angles due to the DM interactions. Dashed 
lines indicate spin directions without angle modulations. 



as 



(2) 



irrespective of the value of J2. 

As a result, the energetical advantage of 6c-cycloidal 
spin state relative to the a6-cycloidal one due to the DM 
interactions, \AE^y^ — A£^g^|, increases as J2 increases. 
The 6c-cycloidal spin state is stabilized when the above 
energy difference dominates over the energetical disad- 
vantage due to the hard-magnetization c axis. The above 
discussion tells us that the a components of DM vectors 
on the out-of-plane bonds are relevant to the stability of 
6c-cycloidal spin state. As seen in Figs. 4(a) and 4(b), 
the diagram for ac=0.38 indeed has a larger regime of 
6c-cycloidal phase than the diagram for ac=0.30. 

The diagram of GTMO is reproduced when ac=0.38, 
and this value is slightly larger than ac=0.30, for which 
the diagram of EYMO is reproduced. This difference 
may be due to different i?-site-radius dependence of the 
GdFeOs-type distortion. According to Ref.,^^ the lattice 
parameters of TbMnOs are equivalent to those of EYMO 
with X ~ 0.85 although TbMnOs is expected to be lo- 
cated at X ~ 0.4 in the T-x diagram in terms of the 
averaged i^-site radius. This indicates that the lattice of 
TbMnOa is more significantly distorted than expected 



from comparison to the EYMO system we may have 

to consider not only the average of R-site radii but also 
their variance. The out-of-plane Mn-O-Mn bond angles 
in GTMO tend to more deviate from 180° than those in 
EYMO, resulting in a larger value of ac- 

There is a slight difference between the experimental 
diagram of GTMO and theoretical diagram for ac=0.38 
compare Fig. 1(b) and Fig. 4(b). In the experimen- 
tal one, the phase boundary between the two cycloidal 
phases slightly bends, and in a narrow region, the system 
exhibits a reentrant behavior with successive transitions 
from be- to ab- and again to 6c-cycloidal phases with low- 
ering T. In addition, the a6-cycloidal spin order is absent 
in the ground state. These points are not reproduced in 
our calculation. This discrepancy may be solved by con- 
sidering effects of /-moments on the R ions, which order 
below ~10 K. By contrast, in the case of EYMO without 
interference from /-moments, the agreement between the 
experiment and the calculation is quite good. 

In i^MnOs, the ab- and 6c-cycloidal spin states are 
stabilized by SIA or DM interaction. Namely, the hard 
magnetization c axis due to SIA gives a relative stability 
to the a6-cycloidal state, while the 6c-cycloidal state is 
stabilized by the DM vectors on the out-of-plane bonds. 
On the other hand, the ac-plane spin cycloid is unfa- 
vorable. Recent first-principles calculations for TbMnOs 
also confirmed this tendency. ^^'^^ 

To summarize, we have studied the spin-cycloidal 
plane flop from ab to be in i^MnOs by constructing a 
microscopic model. We have used realistic parameters 
determined experimentally. We have revealed that the 
flop occurs due to the competition between SIA and DM 
interaction controlled by the GdFeOa-type distortion. It 
has been demonstrated that consideration of this com- 
petition naturally reproduces the experimental T-x di- 
agrams of EYMO and GTMO, and that incorporation 
of SIA and DM interaction is indispensable to describe 
the ME system in i^MnOa. Note that we have also suc- 
cessfully reproduced the sinusoidal collinear spin phase 
in the intermediate T regime, whereas previous theories 
failed. Our finding provides a useful insight for study- 
ing origins and mechanisms of intriguing ME phenomena 
in the MF materials of magnetic-spiral origin. 
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